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I. INTRODUCTION 



The recent results from the v disappearance experiments [l|, |2| have made important 
contributions toward pinning down the elements of the neutrino mixing (PMNS) matrix, 
which we will denote as with elements V^j, a = (e, /i,r), i = (1,2,3). Of the four 
physical parameters in V^, three have now been measured, albeit with different degrees of 
accuracy. Furthermore, some information on the fourth can already be gleaned from the 
known data set obtained by the various extant experiments, as was done in global analyses 
thereof (see, e.g., jl, |^). It seems timely to study the available results in detail, with the 
intent to extract some general properties of Vi, which may be used to suggest directions for 
further investigation. 

In this paper we will concentrate on two aspects of V,^. First, it is interesting to assess the 
impact of the known results on the possible symmetry properties of V^. Next, we address 
the behaviour of \4 as a function of the parameter A = 2\/2G p^eE , the induced neutrino 
mass in matter. This is especially relevant to the long baseline experiments (LBL) (for an 
incomplete list, see, e.g., Ref.js], [6| and the references therein), which are generally regarded 
as the "future" of neutrino physics explorations. 

The analyses of these issues are facilitated by a judicious choice of parameters for V^. To 
begin, in this study we will use the rephasing invariant parametrization introduced earlier. 
It consists of six parameters {xi,yj), which satisfy two constraints. Prior to the recent 
measurement on Ue — ^ t'e, which results in a "large" iKsl, it seems probable that K is — r 
symmetric 0,13], |V^i| = \Vri\. These conditions, when expressed in the variables Qi = Xi+i/i, 
are just Qi = 0. Now that {Vesl is non- vanishing, exact — r symmetry becomes less likely. 
The question remains: " Is there an approximate fi — t symmetry, and how good is it? " . 
It seems natural, then, to interpret f2j as the /x — r symmetry-breaking parameters. As we 
shall see, they satisfy a constraint and thus there are only two independent parameters in 
the set Qi. Also, taken together, the known neutrino data actually constrain all the fij's 
so that none of which deviate much from zero. The other physical parameters we choose 
are Xi = Xi — i/i = \Vei\^ . As we shall see, the set (Xj, f2j) is convenient for studying V^, in 
vacuum as well as in matter.. 

For neutrino propagation in matter, it turns out that Xi, together with Aij = mf — rrip 
satisfy a set of differential equations with respect to A, the induced mass. Experimentally, 
the Xj's and Ajj's are all very well measured in vacuum [A = 0). This means that Xi 
and Ajj are well determined, for all values of A. Separately, Qi satisfy a set of differential 
equations containing Xj and Aij. They can be integrated to obtain Qi{A), with the known 
solution of {Xi, Aij) as inputs. Here, however, the initial values are poorly determined since 
there is only one direct measurement (from atmospheric neutrinos) on Qi — Q2, with large 
errors. Nevertheless, it will be seen that the Qi{A)^s are largely constrained. 

To go forward, our analysis suggests that the most urgent task would be another inde- 
pendent measurement on It will clarify the nature of Vi, to a large extent. At the same 
time, the close correlation between vacuum parameters and those in matter also means that 
LBL experiments with variable A can be extremely useful in the study of V^. 



II. REPHASING INVARIANT PARAMETRIZATION 

As was shown before , one can construct rephasing invariant combinations out of 

elements of a unitary, unimodular (detV^ = +1, so that V* = matrix of cofactors of V), 
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mixing matrix, 

^ijk = ^liV2jV3fc = Rijk — iJ, (1) 

where their common imaginary part can be identified with the Jarlskog invariant J (l6| . 
Their real parts are defined as 

(-R123, -R231, -R312; -R132, -R213, -R321) = {xi,X2, X3; yi, y2, ys). (2) 

These variables are bounded by ±1: — 1 < (xj, < +1, with yj < xi for any They 
satisfy two constraints 

detF = (xi +X2 + X3) - {yi + 1/2 + 1/3) = 1, (3) 

{xiX2 + X2X3 + xgxi) - (1/12/2 + 2/21/3 + ym) = 0. (4) 

In addition, it is found that 

= xia;2X3 - 2/12/22/3- (5) 

Thus, the physical parameters contained in V can be specified by the set (x, y) plus a sign, 



corresponding to J = ±y .P. 

For applications to neutrino physics, it is traditional to label the matrix elements as V^j, 
a = (e, /i,r), i = (1,2,3,). The relations between {x,y) and |V^jp are given by 



W 



IK 



Xl - 2/1 X2- 2/2 2:3 - 2/3 

3^3 - 2/2 a^i - 2/3 a;2 - 2/1 I • (6) 

X2 - 2/3 3^3 - 2/1 a;i - 2/2 



One can readily obtain the parameters (x, y) from by computing its cofactors, which form 
the matrix w with w^W = {detW)I, and is given by 

^ Xl + 2/1 X2 + 2/2 a;3 + 2/3 \ 
w = X3 + 2/2 a^i + 2/3 3^2 + 2/1 • (7) 
\ X2 + 2/3 2:3 + 2/1 Xl + 2/2 / 

Note that the elements of are bounded, — 1 < Wai < +1, and 

'^Wai = Yl = detiy, (8) 

i a 

detw = - Ei/| = E^^ + Eyv (9) 

where the constraint equations Eq. ([3]) and Eq. (jlj) have been used. 

The constraints Eqs. ([3j) and (4) can be easily derived by using the identity ri23r23ir3i2 = 
ri32r2i3r32i. One can obtain other useful relations when we consider product of the form 
^ijk^imn- Thus, the well-known rephasing invariant expression Uff = VaiVpjV*jVpi consists 
of four such terms. For instance, 

n23 = (2/1 - T'J){y2 - iJ) - {x2 - iJ){xz - iJ) 

- {xi - iJ){x2 - iJ) + {x2 - iJ){y3 - iJ)- (10) 
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The combination H^- has the additional property that it is rephasing invariant even if 
det^ = e*^ 7^ 1. Another useful formula (with det\^ = +1) is 

kVaiypiy-yk 

'yk 

= iVajl^lVpil^ + ^eap'yeijkV*jV^iV*^, (11) 

'yk 

where the second term in either expression is one of the F's (r*'s) defined in Eq. (Q. 

We now turn to combinations of the form [yi — iJ){ym — iJ) — {xi — iJ){xj — iJ). As an 
explicit example, consider (?/i — iJ){y2 — iJ) — (a^i — iJ){x2 — iJ) = {yiy2 — X1X2) + iJ{xi + 
X2 — Hi — 2/2)5 or 

VelV,2V:^V^3Vr3 = " 0:1X2) + lJ{l " iKsH- (12) 

In general, for a 7^ /3 7^ 7, i 7^ j 7^ /c, 

VajVakV*iVj3iV^i = {ymUn " XfeXj +iJ{l- \ Vai\^). (13) 

Here, if \Va,i\^ = Xa — yi, then b^c^a, m^n^l. Thus, if we take the matrix elements in 
the a-th row and the i-th column, complex conjugate the vertex (y*^), then the product is 
rephasing invariant and has a well-defined imaginary part. In fact, Eq. f|T3l) provides another 
way to compute J. For instance, in the standard parametrization [l3|, if we take a = e and 
z = 3, we quickly recover the usual expression for J. The real part of Eq. f|T3l) is also useful. 
It enables us to compute other physical variables and, if \Vai\ <C 1, set stringent bounds on 
them. We will discuss these applications in sec. IV and V. 



III. CHOICE OF VARIABLES 

Over the past couple of decades, a wealth of information has been gathered by neutrino 
oscillation experiments. It would be useful to analyze the available data systematically so 
as to gain an overview of the neutrino mixing matrix. To this end it is important to choose 
a set of parameters which can bring out clearly the salient features of K/. In this paper 
we propose to use certain combinations of the variables (x, y) which, as we shall see, can 
highlight the symmetry properties of V^. In addition, they have simple behaviors when used 
in the study of neutrino propagation in matter. 

Specifically, we choose the parameters, 

Xi = Xi-yi = iKiT = Wei, (14) 

= Xi + yi = eijkWf^jWrk = Wei, (15) 
where i = (1, 2, 3). Note that —l<fli<l, and 

W,,-Wr^ = -^e,,k{^,-Qk), (16) 

4J2 = X1X2X3 + Xin2^3 + X2niQ3 + Xafiil^a. (17) 

These variables are considered to be functions of A = 2\/2GFn(,E, the induced neutrino 
mass, which will be used when we discuss neutrino propagation in matter. For the specific 
case of vacuum values, A = 0, we will use the notation and f2°. 
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The variables Xj and fij are not independent. Obviously, 

^X, = l. (18) 
The set Qi also satisfies a simple constraint. From 

WeiWei = detW, (19) 

and, using Eqs. ([3]) and (4), 

detH^ = ^ - ^ y| = 5: X, + ^ y, = ^ 1]„ (20) 

we find 

5]fi,(l-X,) = (21) 

The two constraints Eqs. (fT8|) and (21) are equivalent to Eqs. ([3]) and (4), but Eq. (21) 
is easier to implement since it is linear in Note also that if all the f2j's are equal, as 
happens when Wfj,i = Wri, then Qi = 0. 

Thus, the rephasing invariant parametrization of consists of the set (Xj,fij), subject 
to two constraints, Eqs. (ITSl) and (!2T|) . Together with the mass differences, Aij = Di — Dj, 
Di = , they form a complete set of parameters for the neutrino oscillation phenomenology. 

Before the recent measurements of iKsP, the possibility of a vanishing iKsP and the 
equality W^s = Wt-s led to the hypothesis of /i — r exchange symmetry for neutrino mixing, 
= Wri- With the confirmed small, but non-vanishing, iKsP, A'- — t symmetry becomes 
less likely (although not excluded). Nevertheless, it is interesting to analyze the property of 
Viy under the fi — r exchange operation. To this end, let us introduce a — r parity operator, 
Pfj,r, such that 

Pf.r : W^i ^ Wr^. (22) 

From Xi = 1 - {W^i + W^i) and Eq. ([ISD, we find 

P^r ■■ X, o X„ (23) 
P^, : ^ -a, (24) 

In addition, since a — r exchange does not affect the eigenvalues of the neutrino mass 
matrix, we also have 

P^r ■■ Aij ^ A^j■ (25) 

These results are independent of A, since the matter effect only contributes to the e — e 
element of the effective neutrino Hamiltonian. Finally, from Eq. f lT7|l . we see that is also 
invariant under P^,-, J'^. 

Thus, the neutrino parameters can be classified as 1) even under P^,- : Xj, Ajj, J^; 2) odd 
under P^t- : fij. The quantities fli serve as symmetry-breaking parameters - they provide a 
measure of how good/bad the ^ — t exchange symmetry is. 

We summarize our results in the matrix: 

(XXX 
|[(i-Xi) + (fi3-^^2)] |[(l-X2)+'(^]l-^]3)] \[{i-X;) + {n2-n{)] 
|[(i - Xi) - in, - n,)] |[(i - X2) - (^]l - n,)] |[(i - X3) - {n^ - n,)] 

(26) 

Although in general the parameters fij have the range, —1 < VLi < +1, given the current 
data, it will be seen (Sec. IV and V) that they are all small, \VLi\ ^ O(10~^), while vanishing 
values are not excluded. 
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FIG. 1: At the la level, the measured values for — ^-^e bounded by the dashed lines |3] and 
the dot-dashed lines 0], while the bound on can be represented by the solid lines. It is seen 
that the two parallelograms allowed by Refs. ^ and Q do not overlap at the la level. 
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IV. NEUTRINO MIXING IN VACUUM 



Having settled on the parameter set (Xj, fij), we turn now to the question of their actual 
numerical values. Since the experimental measurements have been given in terms of the 
standard parametrization, we need to transcribe the results into the (Xj, variables. In 
so doing, some informations are bound to be lost in translation. Our numerical results are 
thus only approximate. More precise ones can only be obtained by analyzing directly the 
experiments in terms of the parameters (Xj,f2j). 

For the actual numbers we will use the summaries from existing global analyses 0, 0], 
after making the proper conversion of variables. 

First, the Xj's are all well-determined. There are slight differences between the two 
global analyses. Also, the cases for normal and inverted mass spectra are not significantly 
different. We quote, approximately, X^ = 0.32 ± 0.016 0; X° = 0.31 ± 0.016 fl, and 
X3O = 0.025 ± 0.003 |3|]; X^ = 0.026 ± 0.003 fl. Of course, X? = 1 - X^ - X3O. 

Our knowledge on is far less certain. At the lex level, we have fi^ — = 0-14 to 0.24 
jsj; = —0.14 to 0.08 01 . The discrepancy between these two results, as well as the large 
percentage errors in each, is a reflection of the poor quality of them. 

To complete the list we need one more piece of information on Despite the lack of 
another independent measurement, it turns out that the known values of can already 
set a stringent bound on fig. To see that we return to Eq. f|T2|) in Sec 11, 

VelVeiyt^y^zVrZ = (^1^/2 " 0:1X2) + ^(1 " \Ve3\') . (27) 

This equation is especially useful if iKsP = ^ 1. In this case its LHS is significantly 
bounded since, in general, |l4inv;2|^ < 1/4 and iV^sPlKsP < 1/4. {E.g., 4|l4iP|V;2p = 
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(iKlP + \Ve2?Y - (iKlP - \Ve2?? < !)• NoW, 

2(l/il/2 - = -(Xsfii + Xl^]2). (28) 

Using vacuum values, 

+ X°fi° = (29) 

This follows from I]^i(l ~ Xj) = 0, with the approximation X3 = 0. For V,j in vacuum, 
actually = 2/9 and | V^^gpll/o p ^ I/4. Thus, with = e\ 

and 

< V2e/3 ^ 0.07. (31) 

Also, in the same approximation, 

+ 2fi° + = 0. (32) 

The measured values of Ql — fi^ can now be combined with the above bound ( | i^? + 2^2° | < 
0.21) in a plot in the — Q2 plane, as shown in Fig. 1. Here, the solid line correspond to 
the bound in Eq. fl?T]) . We emphasize that this bound is robust, with possible deviations of 
no more than about 10%, largely from the errors in |V"g3|. On the other hand, the dashed 
and dot-dashed \^ lines have large errors, corresponding to the considerable uncertainties 
in Qi — 0,2- Nevertheless, the allowed regions for and fi^ (and also ^3) are essentially 
confined to the neighborhood of the origin. Despite the significant ambiguities, it seems that 
a fair assessment is given by < O(10~^). 

In summary, given the incomplete results that are available at the present, one can already 
deduce useful and quantitative information about all of the four physical parameters in 
(X°, f2°). Clearly, the next step would be to have a precision measurement on It is most 
useful to concentrate on fig, since this is equivalent to a measurement of J°, according to 
Eq. §0^. 



Experimentally, the determination of fig can be quite challenging. First, it is small (^ 
0.07). Second, its contribution to neutrino oscillation probabilities is not easily disentangled. 
We recall that the probability for z/° — z/^ oscillation is given by 



P(^- ^ jyf^) = S^^ - Re{U:f) sin^ 



± 8(l-5„^)J(sin<l>2isin<l>3isin<l>32), (33) 
where $ij = AijL/AE. The amplitudes Reijlff) are simple functions of Wai, or Xj and fij, 

tQ:/3\ 



according to Eq. flTTl) . and are listed in Table I. In this table, it is sufficient to list i?e(n^- ) 
with {a, 13) = {e,fj.). This is because J2aR^0^?f) = 0? also Re{Illf) = -Re(n^°), since 
U'^f = (n(^")*. For vacuum values, $3^ ~ $32; so we also list separately the combinations 
Re{Ut,^) + Re{U",^). 

In searching for amplitudes that contain fi3, it is clear that, while -Re(n3f ) or Re(Il'^2) 
does depend on fis, once we make the combination Reijl^i + Ilg^), very few amplitudes 
have that property. In fact, the six amplitudes of the form Re{Il2i) and Re{Il'^i + Ilgf ) 
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Reim) WelWe2 

Re{Ull) WelWe3 

i2e(n||) We2We3 

i?e(n|f + n||) Wesil-Wes) 

i?e(n^f) W^iW^2 

i?e(n^f) W^iW^s 

Re{U^t) WeiW^2 - XI 

i?e(n|^2') We2W^3 - X2 

TABLE I: The amplitudes i2e(n^^) are simple functions of Wai, or and ilj 



fall into three groups. 1) ReilV^i) and i?e(n35 + Ilgl): They do not contain fij, and have 
been used to determine successfully. 2) i?e(n3f + 1132) ^^'^ R^^'^i + n32)- Here, 
W,3{l-W,3) ~ i[l-(fi2-^^i)'], and -PVesW^MS - -X3[l + (^2 -^^i)]/2, for X3 « 1. Thus, 
for vacuum values, they can only be used to determine (Q^ 3) The amplitude Reilli^i) 

and Re{Jl2i) do depend on ^3. Substituting in the approximate vacuum values, = 2/3, 
XO ^ 1/3, X3O ^ 0, and dropping quadratic terms in Qi, we find /2e(n^f )° = i-|(2fi°+fi°). 

Thus, to gain access to (actually ^1 + 2^1), one has to isolate the amplitudes Reijli^i)^ 
and i?e(n2i)°, which is by no means easy. We can only hope that the technical difficulties 
involved can be overcome in the near future. 



V. NEUTRINO MIXING IN MATTER 

When neutrinos propagate in matter, their interactions induce a term in the effective 
Hamiltonian, given by {H)ee = A = 2y/2GFneE M. The mass eigenvalues and the mixing 
matrix are now functions of A. It was shown [15| that they satisfy a set of differential 
equations listed in Table II, together with 

These equations are derived from 

{V + dVy{H + dH){V + dV) = D + dD, (35) 

which is just the flavor-basis version of the familiar perturbation theory in quantum me- 
chanics, 

{l + dVy{H + dH){l + dV) = D + dD, (36) 

written in the mass eigenstate basis. Thus, the quantum mechanical result dDi = {i\dH\i) 
becomes Eq. (131|) . dD^ = \Vei\^dA, when dH has only an (ee) element, {e\dH\e) = dA. 
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1/Ai2 I/A23 I/A31 

dA 
dA 

—Q1X3 — ^IsXi + ^2^3 + ^3^2 ^2^3 ~ ^3^2 —^1^3 + ^3^1 



niX2 - ^2X1 -niX2 - ^2X1 + 1^1X3 + nsXi -^1X3 + ^3X1 

^1,1X2 — ^2X1 ^2^3 ~ ^3X2 ^1X2 + ^2X1 — r22^3 — ^^3^2 



TABLE II: The differential equations for Xi and ilj in matter, expressed as the sums of terms 
proportional to 1/Ajj. 



Similarly, the equations for Xi and fij are just rephasing invariant combinations constructed 
out of the well-known formula dVij = {i\dH\j) / [Dj — Di), after its conversion into the flavor 
basis. 

The above arguments should help to illuminate the nature of these differential equations 
and to put on firm grounds the use of them in solving the mixing problem. It is noteworthy 
that the equations for and Xi are independent of These equations are even under 
Pfj,r so that, barring the possible appearance of higher order terms, the f2's are absent. 
Similarly, the differential equations for f2j are odd under P^^r, so that they are all linear in 
Qi. Previously [ll], we studied the A-dependence of Xi and Di, assuming Qi = 0. We now 
see that even with Qi ^ 0, the conclusions remain valid. We only need to add Qi to the list 
and find out how they behave. 

Let us first summarize the results for Xi and Qi. Because the mass differences and 
A32 are widely separated, it is known (see e.g., Ref. 13]) that the three-fiavor problem 
can be well approximated by two, two-fiavor, level crossing problems, occurring at A = Ai, 
when A21 = min, and at A = A^, for A32 = min. In our present formulation this is just the 
pole dominance approximation. Near A ~ Ai, the differential equations for {Di, D2, Xi, X2) 
are reduced to the pole term oc I/-D21, plus Eq. flMl) . Similarly, for A ~ Ah, variations of 
{D2, D3, X2, X3) are governed by Eq. flMl) and the pole term oc I/-D32. Consequently, the 
function Di{A) takes turn to rise linearly, while Xi behaves as step-functions. We list these 
results by dividing A into three regions, 

I) A<Ai: Di oc X^A, D2 oc X^A, D3 = constant, Xi = X\ 



II) intermediate Ai region, A; < A < A^. {Di,D3) = constant, D2 oc A, Xi 
2 1, -^3 — -^3! 



Ill) dense medium Ad, A > Ah. {Di,D2) = constant, D3 oc A, {Xi,X2) 0, X3 ^ 1. 
Detailed graphs for these variables will be given in Sec. VI. 

Having obtained the functions Di{A) and Xi{A), we can now use them to evaluate Qi{A). 
While one can integrate the differential equations for Qi (Table II) numerically, which will 
be presented in Sec. VI, much of the behavior of Qi{A) can already be inferred from the 
constraints (Eq. fl2ip ) and the nature of the differential equations they satisfy. As with 
Xi{A) and Di{A), we will now invoke the pole dominance approximation, so that Qi{A) 
undergo rapid variations only near the resonances A ^ Ai and A Ah. Outside of these, 
Qi{A) are essentially flat. Their values can be summarized as follows. 
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FIG. 2: The parameters X-i (left column) and Xi (right column) as functions of A/Ag^ for both 
the normal (solid) and the inverted (dashed) mass spectra. The initial values in vacuum are well- 
measured, with X? = ~ (2/3)(l - Xl), X^ = X^^ (1/3)(1 - Xl), where X^ = X^ = {V^^l^. 




I) In the intermediate A region, Ai < A < A^. fii -> 0, ^2 — ^ ^2 ~ ^i' ^3 ~^ 0- Here, 
Xi — )■ 0, X2 — )■ 1 — X3 = 1, X3 — )■ X3. Also, a more precise Xi{A) can be obtained by 
using X?X0(A0J2 _ XiX2(A2i)2 (Eq. (4.12) of Ref. O). With X2 ^ 1, (A2i)2 = A^, 
and X°XO = 2/9, we have Xi{A) = {2/9){A^2jAf, for Ai < A < At. Xi{A) is thus already 
very small for A > SAgi. Now we turn to the identity Eq. (|T3|) . for a = e, / = 1, 

V^lVrlV:,Ve2Ve3 = (-X2X3 + ViVs) + ^^(1 - l^^ll'). (37) 

With \Vei\ = e' < 1, ir^iKil < 1/2, = e < 1, and = -(^2X3 + ^3X2) = 

— 2(^2X3 — 112113), we find the bound 

\ni\<ee, (38) 

for Ai A ^ Ah. At the same time, the constraint ~ ^i) = after putting in the 

values of Xj, becomes 

Q-^ + Q.^^ 0. (39) 

Finally, referring to Table II, we have 

±{n2-n^) = o, (40) 
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FIG. 3: The qualitative plots for Di (dashed), D2 (solid), and D3 (dot-dashed) under normal (left) 
and inverted (right) hierarchies. Note that in each plot, the curves in the positive region of A 
represent Di for the z^-sector, while that in the negative region of A represent Di for the z^-sector. 

D, D, 



at the Ai pole. Thus, 

= 1^2-1^1 ^1^2 (^) (41) 

for Ai<A< Ah. 

II) For dense medium, A ^ A^. fii — )■ 0, — ^ 0, — )■ — + ^i- this region, 
Xi — )■ 0, X2 — )■ 0, X3 — )■ 1. Following steps as above, we see first that fl2 0. The 
constraint equation then gives ^1+^2 = 0, and, using the property of the pole at A^, 
d{Q3 — Q2)/dA = 0, we find the results listed. Note that the initial conditions for the 
differential equation are (Aj) = 0, f22(A) = ^2 ~ ^i- 

Our results show that fli{A) behave rather simply as functions of A. The transition 
regions near Ai and Ah, however, are not covered. For these we need to solve the differential 
equations numerically, which will be given in Sec. VI. After these transition regions, the 
resonances act to "quench" many of the parameters. Thus, as A increases, much of the 
information in the original set {X^, Q^) will disappear, and Wu contains fewer and fewer free 
parameters. Physically, this is a consequence of the decoupling of the |e) state, whose mass 
rises roughly as a linear function of A. Starting as a mixture of |1) and |2), it evolves into 
an almost pure |2) state, and eventually becomes just a |3) state. For A ^ m|, the system 
is made up of |e) (|3)) plus a completely decoupled two-flavor system, consisting of and 
|r), or |1) and |2). The vestige of the original /i — t mixing, (fi^ — fi?), however, remains 
undisturbed by the change in the e— channel, coming from Aifge- 

So far we have tacitly assumed the "normal" ordering of the neutrino masses. For the case 
of "inverted" ordering, the resonance at Ah is absent, but the rest of our analyses remain 
valid. 



VI. NUMERICAL SOLUTIONS 

We summarize our analyses of the parameters here with their numerical solutions. As 
indicated earlier, the quantities X^ and A° = Di — Dj are all well measured at A = 0, 
with errors ~ 10%. We show both Xi (for i^) and Xi (for i?) as functions of A/A21 in 
Fig. 2 with the initial values: = ~ (2/3) (1 - ^30) and X^ = X^ ~ (1/3) (1 - X^), 
where X3 = X3 = |Vg3p is taken from both Refs. [3] and It is seen that the paths 
for u and u begin to evolve apart when A ^ Ai. In addition, Xi and X2 are insensitive 
to the mass hierarchy, while each of X2, X3, Xi, and X3 evolves diversely under different 
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FIG. 4: The numerical solutions of 0,i for the i/-sector under both the normal hierarchy (solid) and 
the inverted hierarchy (dashed). Two different sets of initial values are adopted: (05)^2)^3) — 
(0.18,-0.015,-0.05) fl (left column) and (Jl?, Q^, J75]) = (0.06,0.045,0.05) fl (right column). 



ill n, 

0.2 t 0.2 h 




-0.1 - -0.1 



mass hierarchies when A > Ai, where the higher resonance begins to affect. Note that with 
the exchange: normal ■H- inverted, the "trends" of the following curves evolve similarly: 
Xi ~ X2, X2 ~ Xi, and X3 ~ X3. The evolution details of Xi and Xi can be readily 
inferred from how Di varies, as indicated by Eq. f l5^ . As an illustration, we provide the 
well-known qualitative plots of Di in Fig. 3. 

On the other hand, are only roughly known from global fits, with typical errors of 
30 ~ 100%. The evolution equations for f2j in Table II can be expressed in a compact form, 

-— - = ^ — — [5ij{^liXk — i^kXi) — 5ik{i^iXj — ^IjXi)] 

j>k 

+ E -^[-{^^X, + Q,X,) + {Q,X, + QkX,)]. (42) 

The numerical solutions of Qi for both mass spectra are shown in Fig. 4. It is seen that 
the evolution of Qi is sensitive to the choice of initial values, which are not quite settled 
experimentally. The general features, however, can be fairly understood from the constraint, 
Eq. (ETl) . and the nature of the evolution equations, as discussed in Sec. V. Note that the 
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FIG. 5: The numerical solutions of 0,i for the z^-sector under both the normal hierarchy (solid) 
and the inverted hierarchy (dashed). We adopt two different sets of initial values: (i75)^2i^3) = 
(0.18,-0.015,-0.05) 8 (left column) and (0?, J^^, J7[]) = (0.06,0.045,0.05) fl] (right column). 
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curves for the inverted mass ordering follow more closely the trend we discussed in Sec. V 
since they are not distorted by the higher resonance. 

For the z/-sector, the evolution of fli in matter is shown in Fig. 5. The analyses in 
Sec. IV remain valid to the understanding of their general features. Note that Qi (Qi) 
share similar qualitative properties as Xi (Xj), e.g., Qi and ^2 are insensitive to the mass 
hierarchy even when the neutrinos propagate in matter, while the matter effect breaks the 
hierarchy degeneracy for each of the parameters, Q2, ^3, ^1, and Cl^, when A > Ai. 

As one scans through the plots, it is noteworthy that all of the parameters undergo rapid 
changes near the resonance positions, but are otherwise more or less flat. This lends support 
to the use of the pole approximation in Sec. V. The values before and after the transition 
regions also agree with the estimates given in Sec. V. For instance, consider Fig. 4. Here, 
fli drops from = 0.18 to almost zero, when A/A^i varies from about 0.5 to 3. At the 
same time, ^3 changes from -0.05 to ~ 0, while Q2 goes from -0.01 to ~ —0.18, not far 
from ^2 — fi'j' ~ —0.19. From A/Agi ~ 3 to ~ 25, Qi stays nearly constant. The other 
plots can be similarly analyzed. To summarize, unless very high accuracy is demanded, the 
approximation presented in Sec. V should be sufficient for most purposes. 
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VII. CONCLUSION 



The physics of neutrino oscillation is governed by a mixing matrix which is unitary, uni- 
modular, and rephasing invariant. As such (or W^) satisfies a number of self-consistency 
conditions (such as Eq. (|T3|) ) which can help to characterize the matrix. To exploit these 
properties, we propose to use the parameters Xi (= Wei) and Qi (= SijkW^jWrk), which 
satisfy two constraints, Eq. ( 1T8|) and Eq. (12T|) . Physically, the set fli offers a measure of the 
— r asymmetry. These parameters are directly measurable, since the neutrino oscillation 
probabilities P(z/" — )■ z/^) are simple functions of them. This parametrization is summarized 
in Eq. (12^. 

Experimentally, the vacuum values of Xj (X°) are well-determined. But for Qi there 
is only one measurement on (^2° ~ ^1)5 which turns out to be small (O(10~^)), with large 
uncertainties. However, using Eq. f l27|) . one can establish the bound ^2° < 0.07, whose 
validity depends on the approximation X3 ^ 1. There is also a sum rule relating to J°, 
Eq. f l30|) . It is thus urgent to have a precision measurement on Q^, although the task can 
be very challenging. 

Turning to neutrino propagation in matter, we find that the parameters have simple 
dependences on A, the induced mass. To a good approximation, there are two resonance 
regions {A ~ Ai and A ~ A^) where they change rapidly. Outside of these there are three 
regions of A in which all parameters take on values that are nearly constant. These are given 
in detail in Sec. V. Starting from {X'j',Q^) for vacuum {A = 0), the list of free parameters 
gets shorter as A increases. For Ai ^ A ^ Ah, there are only two: ^2(^4) — ^2 ~ ^i) 
Xs^Ai) = Xl For A > Ah, it is down to one: QsiAa) = -(fi^ - Ql). Physically, this 
behavior is a consequence of decoupling, and can be tested in LBL experiments. So far, 
the analysis is done under the assumption of normal hierarchy. For the case of inverted 
hierarchy, similar results are obtained with the omission of the higher resonance. 

In conclusion, given the incomplete knowledge that is now available, it is seen that, with 
the help of consistency conditions and the approximations X^ <^ 1 and ^ A^^, much 
can already be learned quantitatively about Wi,, both in vacuum and in matter. It is hoped 
that our analysis will be helpful toward establishing a comprehensive specification of the 
neutrino mixing matrix. 
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